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A new pressure, volume, temperature, and, 
time (PVTt) primary gas flow standard at 
the National Institute of Standards and 
Technology has an expanded uncertainty 
(k=2) of between 0.02 % and 0.05 %. The 
standard spans the flow range of 1 L/min 
to 2000 L/min using two collection tanks 
and two diverter valve systems. The stan- 
dard measures flow by collecting gas in a 
tank of known volume during a measured 
time interval. We describe the significant 
and novel features of the standard and ana- 
lyze its uncertainty. The gas collection 
tanks have a small diameter and are 
immersed in a uniform, stable, thermostat- 
ted water bath. The collected gas achieves 
thermal equilibrium rapidly and the uncer- 
tainty of the average gas temperature is 
only 7 mK (22 x 10 7). A novel operat- 
ing method leads to essentially zero mass 
change in and very low uncertainty contri- 
butions from the inventory volume. 
Gravimetric and volume expansion tech- 
niques were used to determine the tank 
and the inventory volumes. Gravimetric 
determinations of collection tank volume 
made with nitrogen and argon agree with a 



standard deviation of 16 x 10 VJ^ The 
largest source of uncertainty in the flow 
measurement is drift of the pressure sensor 
over time, which contributes relative stan- 
dard uncertainty of 60 x 1 to the deter- 
minations of the volumes of the collection 
tanks and to the flow measurements. 
Throughout the range 3 L/min to 110 
L/min, flows were measured independent- 
ly using the 34 L and the 677 L collection 
systems, and the two systems agreed with- 
in a relative difference of 150 x 10 . 
Double diversions were used to evaluate 
the 677 L system over a range of 300 
L/min to 1600 L/min, and the relative dif- 
ferences between single and double diver- 
sions were less than 75 x 10 . 
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1. General Description of a PVTt Gas 
Flow Standard 

PVTt systems have been used as primary gas flow 
standards by the National Institute of Standards and 
Technology (NIST) and other laboratories for more 
than 30 years [1 ,2,3,4]. The PVTt systems at NIST con- 
sist of a flow source, valves for diverting the flow, a 
collection tank, a vacuum pump, pressure and tempera- 
ture sensors, and a critical flow venturi (CFV) which 
isolates the meter under test from the pressure varia- 
tions in the downstream piping and tank (see Fig. 1). 
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Fig. 1. Arrangement of equipment in a. PVTt system. 
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The process of making a PVTt flow measurement 
normally entails the following steps: 

1) Close the tank valve, open the bypass valve, and 
establish a stable flow through the CFV. 

2) Evacuate the collection tank volume (Fj) with the 
vacuum pump. 

3) Wait for pressure and temperature conditions in 
the tank to stabilize and then acquire initial values for 
the tank (P^ and Tj). These values will be used to cal- 
culate the initial density and the initial mass of gas in 
the tank (rrij ). 

4) Close the bypass valve and during the "dead-end 
time" when both the bypass and tank valves are fully 
closed, obtain a start time (/'). At the same time, acquire 
the initial pressure and temperature in the inventory 
volume (P^ and 7J'). These values will be used along 
with the equation of state for the gas and the inventory 
volume (Fi) to obtain an initial mass in the inventory 
volume ( m\ ). Shortly after the bypass valve is fully 
closed, open the tank valve. 

5) Wait for the tank to fill to a prescribed upper pres- 
sure (in this system, about 100 kPa), and then close the 
tank valve and obtain the stop time (t^ during the dead- 
end time. At the same time, acquire the pressure and 
temperature in the inventory volume (Pj and T^ ) and 
hence the final mass in the inventory. Open the bypass 
valve. 

6) Wait for stability and then acquire P^ and Tj and 
hence w,,.. 

By writing a mass balance for the control volume 
composed of the inventory and tank volumes (see the 
volume defined by the dashed line in Fig. 1), one can 
derive an equation for the average mass flow during the 
collection time: 



(mj-mjj+(ml-m\j 
m= ; — ; , 



(1) 



t'-e 



or, neglecting the volume changes between the initial 
and final conditions: 



. _v^{p\- p'\)+v,{p[- p\) 



m = 



(2) 



t'-f 



where p is the gas density determined via a real gas 
equation of state: 

PM 

where Z is the compressibility factor, M is the molecu- 
lar weight, and R is the universal gas constant. 

The start and stop times can be chosen at any point 
during the dead-end time as long as the inventory con- 
ditions are measured coincidentally. Why is this true? 



Implicit in the PVTt basis equation [Eq. (2)] are two 
requirements: 1) the measurement of the initial and 
final densities must be coincident with the measure- 
ment of the start and stop times and 2) there must not 
be any other sources or sinks of mass flow to the con- 
trol volume. The second condition is met for the entire 
time that the bypass valve is fully closed, including the 
start and stop dead-end times. It is not necessary that 
the initial and final determinations of the mass in the 
collection tank be done coincidentally with the start and 
stop times because the tank is free of leaks and it is 
advantageous to measure these mass values when the 
tank conditions have reached equilibrium. The freedom 
to choose the start and stop times from within the dead- 
end time intervals allows one to choose times where the 
initial and final inventory densities match, giving 
essentially zero mass change in the inventory volume 
{Am[) and extremely good cancellation of certain corre- 
lated inventory uncertainties. 

Equations (2) and (3) are basis equations for the 
mass flow calculation in a PVTt system and hence are 
the foundation for the propagation of uncertainties 
analysis [5] which will be presented in later sections. 
As will be shown, the most important contributors to 
the uncertainty of the gas flow standard are the collec- 
tion tank volume, the density of the gas in the full col- 
lection tank (both primarily traceable to pressure uncer- 
tainty), and uncertainties related to the measurement of 
the change of mass in the inventory volume. 

We use the volumetric flow unit of L/min throughout 
this publication since the PVTt system is a volumetric 
(as opposed to gravimetric) gas flow standard. The 
L/min flow unit reflects the fact that the PVTt standard 
works over the same range of volumetric flow regard- 
less of the molecular weight of the gas being metered. 
We have used the reference conditions of 293.15 K and 
101.325 kPa when calculating these volumetric flows, 
hence the flow unit corresponds to the widely used 
"standard liter per minute" or slm. 



2. Design and Operation of tlie PVTt 
Standard 

Experience with a previous PVTt flow standard at 
NIST indicated that improvements in temperature and 
pressure instrumentation as well as in the design and 
operation of the new system would be necessary in 
order to achieve the uncertainty goal of 0.05 % or bet- 
ter We now describe aspects of the design and opera- 
tion of the new flow standard that are important to 
achieve this low uncertainty. Section 2.1 describes the 
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measurement of the average temperature of the collect- 
ed gas, Sec. 2.2 describes the procedures that minimize 
the uncertainty of the mass change in the inventory vol- 
ume, and Sec. 2.3 describes the determination of the 
tank and inventory volumes. 



2.1 Average Temperature of the Collected Gas 

One of the most important sources of uncertainty in 
a PVTt flow standard is the measurement of the average 
temperature of the gas in the collection tank, particular- 
ly after filling. The evacuation and filling processes 
lead to cooling and heating of the gas within the vol- 
ume due to flow work and kinetic energy phenomena 
[6]. The magnitude of the effect depends on the flow; 
however, the temperature rise in an adiabatic tank can 
be 10 K or more. Hence, immediately after filling and 
evacuation, significant thermal gradients exist within 
the collected gas. For a large tank, the equilibration 
time for the gas temperature can be many hours. If the 
exterior of the tank has non-isothermal or time varying 
temperature conditions, stratification and non-uniform 
gas temperatures will persist even after many hours. 

In this work, we avoided long equilibration times 
and the difficult problem of measuring the average tem- 
perature of a non-uniform gas by designing the collec- 
tion tanks for rapid equilibration of the collected gas 
and by immersing the tanks in a well-mixed, ther- 
mostatted, water bath (see Fig. 2 and Fig. 3). There are 
two control volumes, a 34 L collection tank and a 677 L 



collection tank. Because the equilibration of the 677 L 
tank is slower, we consider it here. The 677 L tank is 
composed of eight, cylindrical, 2.5 m long, stainless 
steel shells connected in parallel by a manifold. Each 
shell has a wall thickness of / = 0.6 cm and an internal 
radius of a = 10 cm. Because all of the collected gas is 
within 1 cm of a nearly isothermal shell, the gas tem- 
perature quickly equilibrates with that of the bath. 
After the collected gas equilibrates with the bath, the 
gas temperature is determined by comparatively simple 
measurements of the temperature of the recirculating 
water. Remarkably, the water temperature measure- 
ments made with 14 sensors had a standard deviation of 
only 0.4 mK during a typical, 20 minute long, equili- 
bration interval. In Sec. 2.1.1, we describe the bath; in 
Sec. 2.1.2, we discuss the equilibration of the collected 
gas. 

2.1.1 The Water Bath 

The water bath is a rectangular trough 3.3 m long, 
1 m wide, and 1 m high. Metal frames immersed in the 
tank support all the cylindrical shells and a long duct 
formed by four polycarbonate sheets. The duct sur- 
rounds the top, bottom, and sides of the shells: howev- 
er, both ends of the duct are unobstructed. At the 
upstream end of the bath, the water is vigorously stirred 
and its temperature is controlled near the temperature 
of the room (296.5 K) using controlled electrical 
heaters and tubing cooled by externally refrigerated, 
circulated water A propeller pushes the vigorously 
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Fig. 2. A schematic of the PVTt collection tanks, water bath, duct, and temperature control elements. 
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Fig. 3. A photograph of the two PVTt collection tanks submerged in the temperature controlled water bath. 



Stirred water through the duct along the collection 
tanks. When the flowing water reaches the downstream 
(unstirred) end of the trough, it flows to the outsides of 
the duct and returns to the stirred volume through the 
unobstructed, 10 cm thick, water- filled spaces between 
the duct and the sides, the top, and the bottom of the 
rectangular tank. 

The uniformity and stability of the water temperature 
was studied using 14 thermistors. The thermistors were 
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bundled together and zeroed at one location in the 
water bath. Then, they were distributed throughout the 
water bath. Figure 4 plots data recorded at 5 s intervals 
from these 14 thermistors. Nearly all of the data in Fig. 
4 is within ± 1 mK of their mean and the standard devi- 
ation of the data from their mean is only 0.4 mK. 
The largest temperature transients occur where the 
mixed water enters the duct, indicating incomplete 
mixing. The tank walls attenuate these thermal 
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Fig. 4. Temperature data for 14 thermistors distributed in the water bath. 
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transients before reaching the collected gas. Thus, after 
equilibration, the non-uniformity of the water bath and 
the fluctuations of the average gas temperature are less 
than + 1 mK (3 X 10"^ T). 

2.1.2 Equilibration of the Collected Gas 

For design purposes, we estimated the time constant 
(Tgas) that characterizes the equilibration of the gas 
within the collection tank after the filling process. The 
estimate considers heat conduction in an infinitely 
long, isotropic, "solid" cylinder of radius a [7]. For the 
slowest, radially symmetric heat mode, Tg^, = 
(a/2.405)VD-r, where Dj is the thermal difflisivity of the 
gas. This estimate gives Tg^s = 80 s for nitrogen in the 
677 L tank. This estimate for Tg^j is too large insofar as 
it neglects convection, conduction through the ends of 
the tanks, and the faster thermal modes, all of which 
hasten equilibration. The time constants for heat to flow 
from the gas through the tank walls and the time con- 
stant for a hot or cold spot within a wall to decay have 
been calculated and found to be less than a second. 
Therefore, we expect the collected gas to equilibrate 
with a time constant of 80 s or less. 

The equilibration of the collected gas was observed 
experimentally by using the tank as a constant-volume 
gas thermometer. After the tank valve was closed, the 
pressure of the collected gas was monitored, as shown 
in Fig. 5. Our analysis of data such as those in Fig. 5 
leads to the experimental values T^^ of less than 60 s for 
both the 677 L and 34 L tanks, in reasonable agreement 



with the estimates. The measured time constant and 
Fig. 5 show that a wait of 20 minutes guarantees that 
the collected gas is in equilibrium with the bath, within 
the resolution of the measurements. 

The manifold linking the eight cylindrical shells is 
completely immersed in the water bath. Thus, the gas in 
the manifold quickly equilibrates to the bath tempera- 
ture as well. However, each collection system has 
small, unthermostatted, gas filled volumes in the tubes 
that lead from the collection tanks to the diverter 
valves, the pressure transducers, etc. In Sec. 6.2.1, we 
show the possible temperature variations of these 
small, unthermostatted volumes make very small con- 
tributions to the uncertainty of the gas temperature and 
the flow measurements. 



2.2 Mass Change in the Inventory Volume 

2.2.1 Overview and Strategy 

As outlined in Sec. 1, the start time /' and the stop 
time /^used in Eqs. (1) and (2) are chosen to occur dur- 
ing the brief "dead-end times" (< 100 ms) when both 
the tank valve and the bypass valve are closed, i.e., we 
use a "zero overlap" diversion [8]. This choice has the 
advantage of clear mass balance accountability for all 
the gas flowing through the critical flow venturi during 
both diversions and the tank filling. Unfortunately, it is 

f i 

difficult to determine either m^ or w, and hence the 
change in mass within the inventory volume accurately 
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Fig. 5. The equilibration of pressure and temperature immediately following a collection tank filling, 25 L/min 
in the 34 L tank. 
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(especially at high flows) because both the pressure and 
the temperature in the inventory volume rise rapidly as 
the flow through the critical venturi accumulates in the 
inventory volume (see Fig. 6). 
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Fig. 6. Experimentally measured data (25 L/min in the 34 L tank) 
and thermodynamic model predictions for zero and non-zero sensor 
time constants. The model outputs demonstrate that neglect of sensor 
response causes significant error in the measurement of inventory 
conditions. 

Our strategy for dealing with the inventory mass 
change has two elements. First, by design, the invento- 
ry volume Vi is much smaller than the collection tank 
volume Vj. (For the 34 L system, Vj/Vj = 500; for the 
677 L system, VjlVi = 700.) Thus, the uncertainty of 
mass flow is relatively insensitive to uncertainty in m\ 
and m\ since both are small compared with the total 
mass of collected gas. Second, we choose i' near the 
end of the dead-end time and we chose t^ such that 
P{t')=P{t^. These choices define a "mass cancella- 
tion" method: since the initial and final inventory den- 
sities are essentially equal. Am, is nearly zero. In fact, 
we will assume that Awj is zero and consider the quan- 
tity only in terms of flow measurement uncertainty, not 
as part of the flow calculation. Symmetry of the inven- 
tory transients (see Fig. 7) and the mass cancellation 
method also give uncertainty benefits due to high cor- 
relation in the uncertainties of pressure and temperature 
measurements for Awj. 

We tested our strategy for choosing /' and t^ for both 
the 34 L and the 677 L flow standards. (See Sec. 4 for 
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Fig. 7. Superimposed inventory data traces for a start diversion and 
a stop diversion in the 34 L tank at 25 L/min demonstrating "sym- 
metric" diverter valve behavior. The stop dead-end time was approx- 
imately 15 ms longer. 

details of these tests.) To test the 34 L system, we col- 
lected identical flows spanning the range 3 L/min to 
1 00 L/min in both the small and the large tanks, using 
the large tank as a reference for the small tank since its 
inventory uncertainties are quite small in this flow 
range. To test the 677 L collection system, we collected 
identical flows in the 677 L tank following two differ- 
ent protocols. In the first protocol, the inventory vol- 
ume was dead-ended at the beginning and end of the 
collection interval in the usual manner. In the second 
protocol the collection interval was divided into two 
subintervals, which doubled Awj and allowed assess- 
ment of its uncertainty contribution. 

These tests indicate uncertainties due to the invento- 
ry volume that are proportional to flow as would be 
expected based on a thermodynamic model of the 
inventory pressure and temperature transients. If the 
inventory uncertainties are considered to arise from 
uncertainty in the collection time, the inventory mass 
change uncertainty found experimentally for the 34 L 
system was u^ = 4 ms x m (200 x 10~* m for its maxi- 
mum flow). For the 677 L system, single and double 
diversions changed the flow measurement by 75 x 10"^ 
m or less. 

In the remainder of this section, we describe condi- 
tions within the inventory volume during the dead-end 
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times using both a model and measurements. The meas- 
urements show that T{i) and P{i) are nearly the same 
during the start and stop dead-end times. Finally, we 
show that Awj is insensitive to the exact choice of t\ 
provided that the condition P{t^) = P{t^ is applied near 
the end of the dead-end time. 

2.2.2 Conditions Within the Inventory Volume 

Figure 6 displays the time dependent temperature 
T{t) and pressure P{t) in the inventory volume of the 
smaller collection system at a typical collection rate 
{m = 25 L/min; collection time = 82 s). The triangles 
(t = 0) in Fig. 6 were calculated from the lumped- 
parameter, thermodynamic model developed by Wright 
and Johnson [6]. The model assumes a constant mass 
flow m at the entrance to the inventory volume. The 
model neglects heat transport from the gas to the sur- 
rounding structure and non-uniform conditions, such as 
the jet entering the volume. For Fig. 6, T{t) and P{t) 
were calculated on the assumption that the diverter 
valve reduced the flow linearly (in time) to zero during 
the interval -0.02 s < t < 0. Experimentally measured 
values of T(t) and P(t) recorded at 3000 Hz (smooth 
curves) are also shown in Fig. 6. Most of the differ- 
ences between the measured curve and the (t = 0) cal- 
culated triangles result from the time constants of the 
sensors used to measure T{t) and P{t). This is demon- 
strated by the agreement between the experimental 
curve and the model results when time constants are 
incorporated (circles). 

In Fig. 6, the calculated curves do not display fea- 
tures that mark either the onset or the completion of the 
diverter valve closing. Thus, even T{t) and P{t) data 
from perfect sensors cannot be used to mark these 
events. For this reason, the times t^ and ;^ were chosen 
at times that were clearly within the dead-end time 
intervals. 

Figure 6 shows that the measured values of T{t) and 
P{f) are consistent with the Wright- Johnson model for 
the inventory volume after allowance is made for the 
response times of the sensors. The consistency shows 
that the behavior of the inventory volume is understood 
semi-quantitatively. However, this is not sufficient to 
accurately calculate the density p{t) from measure- 
ments of T{t) and P{t) because the fraction of the flow 
collected as the valves are closing cannot be deduced 
from the measurements. Instead, we relied on the pres- 
sure sensor to choose t\ The pressure sensor is pre- 
ferred to the temperature sensor because it responds 
more quickly and also because it responds to the aver- 
age conditions throughout the inventory volume rather 



than the conditions at only one location. We choose i' 
near the end of the dead-end time, where the P{t) meas- 
urements are nearly parallel to the T = model. In this 
regime, the derivative dP/dt is large and its dependence 
on precisely how the valve closed is small. Because the 
dependence on how the valve closed has decayed, we 
expect that P{t) will be the same during the start and the 
stop dead-end times, improving the mass cancellation 
as well as the correlation of initial and final inventory 
density uncertainties. 

2.2.3 Near Symmetry of Start and Stop Behavior 
of P(t) 

Figure 7 shows records of T{t) and P{t) taken during 
the dead-end time intervals at the start and the stop of a 
single flow measurement. As before, the data were 
recorded at 3000 Hz for 500 ms and the plots were dis- 
placed along the horizontal axis until they nearly over- 
lapped. The pressure and the temperature at the begin- 
ning of the start dead-end time were slightly lower than 
those at the stop dead end time; however, the two 
records match closely during the dead-end time. This 
implies that the time-dependent densities p(t) also near- 
ly match. 

At both diversions shown in Fig. 7, valve trigger sig- 
nals were gathered along with the temperature and 
pressure measurements using a commercially manufac- 
tured data acquisition card (see Fig. 8). The trigger sig- 
nals originate from an LED/photodiode pair and a flag 
on the valve actuator positioned so that the circuit out- 
put rises to a positive voltage when the valve is closed. 
These valve signals are used to trigger timers which 
give the approximate collection time. 

As represented in Fig. 8, the inventory record is post- 
processed by the controlling program to obtain both the 
initial and final measurements of pressure and temper- 
ature in the inventory volume as well as the final col- 
lection time. A "match pressure" [P{t^)] is chosen that 
falls late within the dead-end time. The match pressure 
value is found in both the start and stop data series and 
the time differences between the match pressure meas- 
urement and the start and stop trigger signals (A;' and 
At^ are determined from the data record. The appropri- 
ate time correction is added to the approximate collec- 
tion time from the timers. Thus, by adjusting the collec- 
tion time using the inventory data records, the initial 
and final inventory pressures and temperatures are 
nearly matched, leading to nearly equal initial and final 
inventory densities and inventory mass cancellation. 
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Fig. 8. Data records of inventory sensors and valve trigger signals are used to adjust the collection time and improve cancellation of 
the initial and final inventory mass as well as inventory uncertainties. 



2.2.4 Insensitivity of Am, to the Match Pressure 

Figure 9 shows the total correction time as a function 
of the match pressure for two flows in the 34 L system. 
The 100 L/min flow is very high for the 34 L tank, hav- 
ing only an 18 s collection time. Match pressure is 
shown as a percentage of the range of pressures meas- 
ured during the diversion transient. For a perfectly fast 
system (valves and sensors), these plots would be hori- 
zontal lines, i.e., any chosen match pressure would 
result in the same time correction. However, for the real 
system with its inevitable limitations, the match pres- 
sure does matter Exploring the possible reasons for this 
is valuable for improving the system and for obtaining 
an accurate uncertainty analysis. 
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Fig. 9. The collection time correction vs the match pressure used in 
the inventory mass cancellation algorithm. 

First recall that the inventory sensors have non-zero 
time constants and therefore the measurements they 
provide are damped versions of the real conditions and 
further, the values they report at any given instant are 
subject to the recent history of the pressure or tempera- 



ture value. Second, realize that perfect symmetry of 
conditions before and during diversion is unobtainable 
and that these imperfections and the significance of the 
sensor damping increase with the flow. For example, at 
high flows, the rate of change of pressure during the 
tank filling process is large and it becomes more diffi- 
cult to make the pressure at which the stop diversion 
begins closely match the pressure at which the start 
diversion began (due to sensor response and valve con- 
trol delays). This "trigger pressure difference" will be 
considered again in Sec. 4. As another example, the 
bypass and tank valves may not close at the same 
speed. 

Analysis of the thermodynamic model of the inven- 
tory and its sensors shows that times later in the dead- 
end time give better mass cancellation under these cir- 
cumstances since the sensor output enters a period with 
nearly constant slope that is equal to the real pressure 
slope. The experimental results given in Fig. 9 support 
this assertion: match pressures between 50 % and 90 % 
result in nearly constant correction times, while low 
match pressures (early in the dead-end time) give much 
larger corrections. Based on this analysis, a match pres- 
sure of 80 % has been selected for use in the flow stan- 
dard. Figure 9 demonstrates the insensitivity of Am, to 
a wide range of match pressure values. 

Figure 9 also illustrates the concept that uncertainties 
related to the inventory volume can be treated not only 
as mass measurement uncertainties, but as time meas- 
urement uncertainties as well. One can consider the 
uncertainty in the measurement of time between condi- 
tions of perfect mass cancellation, or one can consider 
the uncertainty in the measurement of inventory mass 
differences between the start and stop times. Both per- 
spectives offer insight and verification of the uncertain- 
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ties of the inventory volume and flow diversion 
process. 



2.3 Measurement of the Tank and Inventory 
Volumes 

2.3.1 Gas Gravimetric Method 

The volume of the 677 L tank was determined by a 
gas gravimetric method. In this method, the mass of an 
aluminum high pressure cylinder was measured before 
and after discharging its gas into the evacuated collec- 
tion tank. The change in mass of the high pressure 
cylinder and the change in density of the gas in the col- 
lection tank were used to calculate the collection tank 
volume. Nominally, 



V = 

grav 



w„ -m, 



Pi- 



'' — V 



(4) 



where the m^ indicates the mass of the high pressure 
cylinder and Vg^^ represents the extra volume tem- 
porarily connected to the tank for the purpose of intro- 
ducing the gas from the cylinder to the tank (usually a 
small volume of tubing and a valve body). The extra 
volume is calculated from dimensional measurements 
or measured by liquid volume transfer methods. 

In practice, a more complex formula than Eq. (4) was 
used to account for a small amount of gas that enters 
the control volume from the room when the cylinder is 
disconnected from the collection tank since the final 
tank pressure was less than atmospheric. For the vol- 
ume determinations performed for the 677 L tank, the 
effect amounts to only 5 x 10"^ Vj. 

The volume determination was conducted with both 
nitrogen and argon gas. In both cases high purity gas 
was used (99.999 %) and care was taken to evacuate 
and purge the system to minimize composition uncer- 
tainties. When nitrogen was used, the aluminum cylin- 
der weighed approximately 4200 g when filled at 12.5, 
and approximately 3800 g after it was emptied to 
55 kPa. When argon was used, the initial and final 
masses were 4440 g and 3820 g, respectively. The stan- 
dard deviation of the six volume measurements (four 
with nitrogen, two with argon) was 1 6 X 10"^ Vj. 

The initial and final masses of the gas cylinder were 
measured using a substitution process with reference 
masses and a mass comparator enclosed in a wind 
screening box. The comparator has a frill scale of 10 kg 
and resolution of 1 mg. The cylinder and a set of refer- 



ence masses of nearly the same weight were alternated 
on the scale five times. The zero corrected scale read- 
ings were then calibrated to the reference masses and 
buoyancy corrected via the following formula: 



w„ 



-m„ 



I — 



~^ Pair 'ext ' 



(5) 



where S represents the scale reading, ref indicates the 
reference masses, p^ is the ambient air density where 
the measurements were conducted, and V^t is the exter- 
nal volume of the high pressure cylinder and its valve 
and fittings. The density of the ambient air was calcu- 
lated from the barometric pressure, the temperature and 
humidity inside the wind screen, and an air density for- 
mula that includes humidity [9]. 

The external volume of the high pressure cylinder 
appears in Eq. (5) due to air buoyancy corrections. The 
external volume of the cylinder was measured by 
Archimedes principle, i.e. by measuring the change in 
apparent mass of the object in two media with differing 
and known densities. One of the media was distilled 
water, and the cylinder apparent mass in the water was 
measured as follows. Liquid was added to the cylinder 
interior until it was nearly neutrally buoyant in the tank 
of distilled water The addition of liquid inside the 
cylinder has no effect on its external volume. The tem- 
perature of the distilled water was raised or lowered 
(thereby changing the density of the distilled water) 
until the cylinder was essentially neutrally buoyant. At 
this point, the apparent mass in the distilled water is 
zero. The temperature of the distilled water was record- 
ed and its density was calculated via an equation from 
the literature [10]. Hence, the temperature of the dis- 
tilled water was used in place of a weigh scale to meas- 
ure the apparent mass in water The apparent mass of 
the cylinder in air (with the liquid still inside) was 
measured using the comparator described above. The 
density of air with humidity was calculated as previous- 
ly described. The external volume of the cylinder was 
calculated for the nominal room temperature (T,^^ of 
296.5 K with the following formula: 



'ext V^ref / 



A A 

m,:, - w,„„, 



[Pwa.Jl+3a(7;,er " 7:ef )] " Pa. [1 +3«( ^ " Tif )] ] ' 

(6) 

where the superscript A indicates apparent mass and a 
is the coefficient of thermal expansion for the alu- 
minum tank. The terms containing a correct for 
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changes in the cylinder volume due to differences 
between the water temperature, the air temperature, and 
the reference temperature. However, for this particular 
case, these thermal expansion issues could have been 
neglected since both the water temperature and air tem- 
perature never differed from T^fhy more than 1.5 K. 
The thermal expansion corrections to the external vol- 
ume were less than 0.5 cm^ or 100 x 10"* V^^ and the 
external volume has a small sensitivity coefficient in 
the collection tank volume determination process. 

The expansion of the external cylinder volume as a 
function of its internal pressure was not negligible. The 
Archimedes principle measurements showed a volume 
increase from 4697.5 cm^ to 4709 cm^ between the 100 
and 12.5 Mpa pressures. This change agreed well with 
predictions based on material properties, and the appro- 
priate experimental values for external volume were 
used in the cylinder mass calculations [Eq. (5)], 
depending on whether the cylinder was empty or full. If 
this issue were neglected, it would lead to relative 
errors in the mass change measurements of about 35 x 
10^. 

2.3.2 Volume Expansion Method 

The 34 L collection tank volume, the inventory vol- 
ume for the large collection tank, and the small inven- 
tory volume were all determined with a volume expan- 
sion method. In this method, a known volume is pres- 
surized, the unknown volume is evacuated, a valve is 
opened between the two volumes, and the density 
changes within the two volumes are used to calculate 
the unknown volume. Applying conservation of mass 
to the system of the two tanks yields: 



ip2-pl) 



V 

extra ' 



(7) 



where the subscripts 1 and 2 refer to the known and 
unknown volumes, respectively. As before, the density 
values are based on pressure and temperature measure- 
ments of the gas within the volumes and gas purity 
issues must be considered. Note that in many cases the 
final densities can be considered the same in both vol- 
umes 1 and 2, but for the determination of the 34 L tank 
volume, elevation differences between the two tanks 
required a head correction to the pressure measure- 
ments and therefore the two densities were not strictly 
equal. The difference in elevation resulted in a relative 
difference in gas density of 20 x 10~* even though the 
two tanks were connected. 



3. Mass Flow Uncertainty 

The uncertainties contributing to the mass flow 
measurement have been quantified in Appendix A and 
now they will be combined. The uncertainty for flows 
between 20 L/min and 2000 L/min of nitrogen or argon 
in the 677 L tank is given in Table 1. The standard 
uncertainty of each sub-component is given in both rel- 
ative (x 10"*) and dimensional forms. The units of the 
dimensional values are given in the third column. The 
relative contribution of each sub-component to the 
combined uncertainty is listed in the fourth column. 
This contribution is the percentage of the squared indi- 
vidual component relative to the sum of the squares of 
all sub-components. The uncertainty from the invento- 
ry volume, the combined uncertainty, the expanded 
uncertainty, and the uncertainty contributions are given 
as a range covering the minimum to maximum flow. To 
calculate their relative uncertainty in Table 1, the tank 
initial density was normalized by the tank final density 
and the inventory mass change was normalized by the 
total mass collected. 



Table 1. Uncertainty of nitrogen flow measurement with the 677 L standard 



Uncertainty category 






Standard uncertainty 


Contrib 


Comments 










(i = 


= 1) 








Flow (677 L, N2) 




Relative (x 


10'^) 






Absolute 


(%) 




Tank volume 




71 








48.44 cm^ 


50 to 23 


see Table A12 


Tank initial density 




10 








2.27xlO"'Vcm^ 


ItoO 




Tank final density 




68 








7.82xl0"^g/cm^ 


45 to 21 


see Table A6 


Inventory mass change 




to 109 








0.084 g 


0to53 


see Table Al 8 


Collection time 




5 








0.3 ms 


Oto 1 


see Table A8 


Std deviation of repeated 


meas. 


20 








0.001 g/s 


4 to 2 




RSS (combined uncertamty) 


102 to 150 














Expanded uncertainty {k - 


= 2) 


204 to 300 
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At the highest flow, uncertainty contributions are 
principally divided between the tank volume, the final 
gas density, and the inventory uncertainty. The uncer- 
tainty falls to 200 X 10"^ m for the smallest flows as the 
uncertainty contributions of the inventory volume 
become negligible. For an air flow measurement, the 
uncertainty of the 677 L system is less than 500 x 10"^ 
m over the entire flow range and the uncertainty is driv- 
en by the tank final density measurement (80 % con- 
tributor). 

Table 2 presents the uncertainty of flow measure- 
ments from the 34 L system for flows between 1 L/min 
and 100 L/min. The expanded uncertainty varies 
between 270x10^ m and 440 x 10^ w. At high 
flows, the significant uncertainty sources are the tank 
volume, the tank final density, and the uncorrelated 
inventory uncertainties. For low flows, the major con- 
tributors are tank volume and final gas density. For air 
flow measurements, the 34 L system has a nearly con- 
stant uncertainty over its entire flow range of about 
500 X 10"^ OT and it is driven by the uncertainty of the 
final gas density. 



4. Experimental Validation of the 
Uncertainty Analysis 

The performance of the flow standards and the valid- 
ity of their uncertainty analyses (and particularly the 
assumed values for the uncorrelated inventory uncer- 
tainties) were checked by performing two types of 
experiments described in Sees. 4.1 and 4.2. 



4.1 Comparison of the 34 L and 677 L Flow 
Standards 

To test both systems, we performed comparisons 
between the two flow standards over the range of flows 

Table 2. Uncertainty of nitrogen flow measurement with the 34 L standard 



where they could both be used: 3 L/min < w < 110 
L/min (3.6 g/min to 138 g/min of nitrogen). The collec- 
tions ranged from as short as 1 8 s to more than 4 h. 

Figure 10 shows the difference in the discharge coef- 
ficients of several critical flow Venturis as measured by 
both the 34 L and 677 L systems, plotted versus flow. 
The agreement between the two flow standards is 
150 X 10"^ m over the entire range tested. The throat 
diameters of the Venturis used for the comparisons 
ranged between 0.3 mm and 1.7 mm. The comparisons 
were done with the same pressure and temperature sen- 
sors associated with the venturi during the testing on 
both flow standards in order to reduce some possible 
sources of discharge coefficient differences. Numerous 
collections were made for each tank at each flow to 
confirm stability of the conditions at the meter under 
test. 

How well should the two systems agree? The differ- 
ence between the discharge coefficients measured by 
the two PVTt systems should be less than the root sum 
square (RSS) of the uncertainties of the two standards, 
especially when one considers that the uncertainties 
due to pressure and temperature measurements are cor- 
related between the two standards. For the lowest flows 
of the comparison range, the uncertainties originating 
from the inventory volume are quite small for both sys- 
tems and the observed differences between them are 
dominated by tank volume uncertainties. From Fig. 10 
it can be seen that the two systems differ by about 
100 X 10^ OT for flows less than 20 L/min. The RSS of 
the two relative standard volume uncertainties is 
137xl0^(;c=l). 

At the higher flows of the comparison range, the 
uncertainties associated with the transient conditions in 
the inventory volume should be negligible in the 677 L 
system, but growing with increasing flow for the 34 L 
system. Because the collection times were l/20th as 
long when using the smaller tank, any timing error (or, 
equivalently any imperfection of the mass cancellation 



Uncertainty category 






Standard uncertainty 


Contrib 


Comments 










ik = 


= 1) 








Flow (34 L, Nj) 




Relative (xlO") 






Absolute 


(%) 




Tank volume 




116 








3.95 cm^ 


72 to 28 


see Table A 13 


Tank initial density 




10 








2.27 X 10"'^ g/cm^ 


ItoO 




Tank final density 




68 








7.82xl0"^g/cm^ 


25 to 10 


see Table A6 


Inventory mass change 




to 170 








0.007 g 


0to61 


see Table Al 9 


Collection time 




15 








0.3 ms 


OtoO 


see Table A8 


Std deviation of repeated 


meas. 


20 








4 X 10"^ g/s 


2tol 




RSS (combined uncertainty) 


137 to 219 














Expanded uncertainty (k = 


= 2) 


274 to 438 
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Fig. 10. Difference in the discharge coefficient of critical flow Venturis calibrated on both the 34 L 
and 677 L flow standards versus flow and the inverse of the collection time for the 34 L tanlc. Also 
plotted is a linear best fit of the data. 



technique) was 20 times more important when using 
the smaller tank. Figure 10 suggests that the inventory 
uncertainties cause the flow standard to read too high as 
the flow is increased, changing the difference by about 
200 X 10"^ m over the range of flows compared. This 
value is comparable to the relative uncertainty of 135 x 
10"^ (A: = 1) contributed by the inventory volume at the 
highest flows in the 34 L system (61 % of 219 X 10^ 
from Table 2). Therefore, the differences observed in 
the comparison are in reasonable accordance with the 
uncertainty analysis. 

Figure 10 also shows the tank comparison results 
from the perspective of time measurement uncertainty 
rather than the mass. We interpreted the comparison 
results using the simplified model m = mit, where m is 
the mass collected and t is the collection time using the 
34 L tank. Making the assumption that there exists a 
constant error in the mass measurement, 5m, for all 
flows (say due to a tank volume error) and a constant 
error in time measurement, 5t, for all flows (say due to 
the time constant of the inventory pressure sensor), we 
can derive a simple linear model for the error in mass 
flow measurement, i.e.. 



5m _ 5m 5t 

m m t 



(8) 



Since the mass collected and the mass error are both 
essentially constant, this model suggests that a linear 
function of the inverse of the collection time should fit 
the tank comparison data. Such a linear function fits 
reasonably well, as shown in Fig. 10, and its slope 
implies a constant timing error of 4 ms. A timing error 
of this order is not surprising because the timing is 
based on a pressure sensor with a time constant of 
approximately 20 ms. The standard deviation of a 
measurement from the best fit line is 24 X 10"^ m. 

The comparison results could be the basis for correc- 
tions to the 34 L system. The intercept of Fig. 10 could 
be used to change the volume of the 34 L tank and 
improve the agreement between the two systems at low 
flows. The slope could be used to improve agreement at 
higher flows. This approach offers the possibility of 
reducing the comparison differences to zero with a 
standard deviation of 24 x 10"^ m. However, these cor- 
rections have not been made at the present time for sev- 
eral reasons. The comparison results are consistent with 
the uncertainty analysis. Also, despite the success of 
Eq. (8), we feel that the inventory uncertainties are 
more related to pressure and temperature than time, so 
it is more appropriate to make improvements in those 
measurements to reduce the slope in the comparison 
data. The volume (or offset) differences can be 
improved by repeating and refining the volume expan- 
sion process used to determine the 34 L tank size. This 
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approach adheres to the definition of primary standard 
for both systems since neither one has been cahbrated 
by a flow measurement against some other flow stan- 
dard. 

During some of the comparison flows, we noticed 
that the pressure downstream of the critical flow ven- 
turi was significantly higher in the 34 L system than in 
the 677 L system (108 kPa vs 100 kPa) due to the 
smaller tube size and resultant higher pressure drop. 
For some of our Venturis (with relatively short dif- 
fusers) this pressure difference caused slight changes in 
the upstream pressure (and the discharge coefficient), 
even at conditions well above the critical pressure ratio. 
Therefore, our assumption that for the same throat 
Reynolds number, the discharge coefficient of the ven- 
turi is independent of the downstream pressure may not 
be perfectly valid. We suspect that some of the differ- 
ences observed between the tanks in Fig. 10 are due to 
the Venturis (even though long diffusers were used). 

In one series of experiments, the trigger pressure dif- 
ference was purposely varied over a range from -2 kPa 
to 27 kPa at a constant flow of 82 L/min in the 34 L 
system. The purpose of the test was to measure the 
dependence of the venturi discharge coefficient on the 
trigger pressure difference and hence assess its influ- 
ence on the inventory volume uncertainties. The tests 
showed a relative change of 10x10"^ in discharge 
coefficient for each 1 kPa change in the trigger pressure 
difference. Since the largest trigger pressure difference 
is less than 3 kPa in the present system, this effect is 
expected to contribute only 30 X 10~* rh to the flow 
uncertainty. Therefore the major contributor to the 
inventory uncertainty appears to be spatial inconsisten- 
cy of the pressure and temperature fields between the 
start and stop diversions or some other, unknown flow 
dependent uncertainty source. 



4.2 Multiple Diversions in the 677 L Flow 
Standard 

To confirm that the uncertainty analysis for the 
inventory volume of the 677 L collection system was 
reasonable, we performed CFV calibrations at identical 
flows following two different protocols. In the first pro- 
tocol, the inventory volume was dead-ended at the 
beginning and end of the collection interval, in the 
usual manner In the second protocol, the collection 
interval was divided into two subintervals, i.e., each 
flow measurement had two start and stop diversions. 
The intermediate dead-end times were set up so that the 
pressure transients in the inventory volume still permit- 



ted the mass cancellation procedure. Breaking the col- 
lection into two subintervals has the effect of doubling 
the uncertainty contribution from the inventory vol- 
ume. The CFV discharge coefficients from the two pro- 
tocols were compared to assess the magnitude of the 
uncertainties introduced by the inventory volume and 
the flow diversion process. Three flows between 300 
L/min and 1 600 L/min were tested and the differences 
in discharge coefficient were all less than 75 x 10~* m as 
shown in Table 3. 

Table 3. Differences in CFV discharge coefficients {C^ for two and 
one diversion in the 677 L flow standard 



Flow (L/min) 

300 

700 

1600 



[Q (2 diversions) - CJ/Q x lO" 

53 ±25 
-27 ±31 
75 ± 122 



5. Conclusions 

The design of a new gas flow standard composed of 
two PVTt collection tanks (34 L and 677 L) has been 
presented. The system is designed to calibrate critical 
flow Venturis for flows from 1 L/min to 2000 L/min. 

The flow standard has several novel features. The 
collection tanks are immersed in a water bath that 
matches the nominal room temperature and is stable 
and uniform to 1 mK. The collection tanks are divided 
into sections of small enough diameter that the gas 
inside them achieves thermal equilibrium with the sur- 
rounding water bath in 20 min or less. This reduces the 
contribution of temperature to the flow measurement 
uncertainty to a low level. 

Uncertainties related to the inventory volume and the 
diversion of gas into the collection tank at the start and 
stop of a flow measurement have been studied in great 
detail. A thermodynamic model of the inventory vol- 
ume during diversion has been utilized to understand 
the large pressure and temperature transients and the 
importance of sensor time constants on the flow meas- 
urement uncertainty. The flow standard is operated to 
achieve "mass cancellation" in the inventory volume, 
thereby taking advantage of correlated sensor uncer- 
tainties to minimize uncertainty contributions from the 
inventory volume. The uncertainty contributions of the 
inventory volume have been considered from both the 
time and mass perspectives. 

The volumes of the collection tanks were measured 
by two methods, a gas gravimetric method and a vol- 
ume expansion method. Six gravimetric determinations 



33 



Volume 108, Number 1, January-Febraary 2003 

Journal of Research of the National Institute of Standards and Technology 



of the 677 L collection tank volume made with nitrogen 
and argon agree with a standard deviation of 16 x 10"^ 

A detailed uncertainty analysis for the gas flow stan- 
dard has been presented in Appendix A. The analysis 
started with the basis equation utilized to calculate flow 
in a PVTt system and followed the propagation of 
uncertainties method suggested by international stan- 
dards. The uncertainties of the sub-components have 
been examined at a fundamental level. 

The uncertainty analysis shows that the 677 L system 
measures mass flow with an uncertainty between 200 x 
10"^ lit and 300 X 10"^ m for a pure gas like nitrogen or 
argon. The higher uncertainty applies to higher flows as 
the inventory transients and the related uncertainties 
grow larger. For the 34 L tank and pure gases, the 
uncertainties range from 440 X 10~^ m to 270 X 10"^ m. 
The uncertainties are larger for the 34 L tank because 
the tank volume uncertainty is relatively greater and the 
inventory volume is relatively larger for the small sys- 
tem. For pure gas measurements, the largest sources of 
uncertainty can be traced to pressure measurement 
(about 70 X 10"^ P) which is the major contributor to 
gas density and tank volume uncertainties. For air flow 
measurements using gas from the existing compressor 
and drier, mass flow uncertainties are about 500 x 10"^ 
rh for both standards and the major contribution is the 
uncertainty in the moisture content of the air 

Comparisons between the 34 L and 677 L standards 
from 3 L/min to 100 L/min show agreement within 
150 X 10'^ m or better Experiments using single diver- 
sions (normal operation) and double diversions to the 
collection tank were used to validate the uncertainty 
estimates of the 677 L inventory volume and the differ- 
ences between these two methods were less than 
75 X 10"^ m. The evaluation results along with compar- 
isons to previously existing gas flow standards support 
the uncertainty statements for the new standards. 

There are opportunities for improvement in the 
uncertainty of the new PVTt flow standards. Pressure 
uncertainties are the most significant contributors, 
through the tank volume determinations as well as 
through the final gas density measurement. Therefore, 
tank pressure sensors with better calibration stability 
over time will be sought. Faster inventory pressure sen- 
sors would improve the mass cancellation procedure 
and reduce the inventory uncertainties. 



6. Appendix A. Uncertainty Analysis 

The uncertainty of a mass flow measurement made 
with the PVTt standard will now be considered follow- 
ing the propagation of uncertainties techniques 
described in the ISO Guide to the Expression of 
Uncertainty in Measurement [5]. The process identifies 
the equations involved in the flow measurement so that 
the sensitivity of the final result to uncertainties in the 
input quantities can be evaluated. The uncertainty of 
each of the input quantities is determined, weighted by 
its sensitivity, and combined with the other uncertainty 
components to arrive at a combined uncertainty. 

As described in the references [5,11], consider a 
process that has an output, y, based on N input quanti- 
ties, Xj. For the generic basis equation: 



y ~ yv-^i^ ■^2'""' -^n)' 



(Al) 



if all the uncertainty components are uncorrected, the 
standard uncertainties are combined by root-sum- 
square (RSS): 



«.w=JS 



JLf;).A 



dy 



3x 



u (x,. ), 



(A2) 



where m(x,) is the standard uncertainty for each of the 
inputs, and u^(y) is the combined standard uncertainty 
of the measurand. The partial derivatives in Eq. (A2) 
represent the sensitivity of the measurand to the uncer- 
tainty of each input quantity 

In cases where correlated uncertainties are signifi- 
cant (as in the following analysis), the following 
expression should be used instead of Eq. (A2): 

uAy) = 



n^j 



dx dx. 



(A3) 



where r(x„ x,) is the correlation coefficient, ranging 
from -1 to 1, and equaling zero if the two components 
are uncorrected. As will be seen in the following analy- 
sis, some uncertainty components in the present system 
are correlated and this leads to a significant improve- 
ment in the uncertainty of the measurand. 
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A simple example of a correlated uncertainty is illus- 
trative. Suppose that a thermometer was used to meas- 
ure a temperature difference. Also suppose that the only 
uncertainty in the thermometer measurement was an 
unknown offset in its calibration. When the difference 
between two temperatures was calculated from meas- 
urements made with this thermometer, the offset would 
cancel and would not contribute to the uncertainty of 
the temperature difference. In this case, the subtraction 
process used to calculate the difference leads to sensi- 
tivity coefficients of opposite sign for the two tempera- 
ture measurements. Since the sensor always has the 
same offset, the uncertainties are perfectly correlated 
[r{Xj, Xj)=\]. When this hypothetical scenario is 
processed through Eq. (A3), the uncertainty of the tem- 
perature difference is zero. Of course in a real case, 
there would be other, uncorrected uncertainties that 
would make the uncertainty of the temperature differ- 
ence non-zero. Nonetheless, the example demonstrates 
that under certain circumstances, correlated uncertain- 
ties will reduce the uncertainty of a measured quantity. 

In the following uncertainty analysis, correlated and 
uncorrected uncertainties will be treated as separate 
components, even if they are related to the same phys- 
ical quantity. For instance, there will be a correlated as 
well as an uncorrected inventory pressure component. 
In this manner, the correlated components can be con- 
sidered as having a correlation coefficient of 1, while 
uncorrected components have correlation coefficients 



of 0. This approach simplifies the process to deciding 
which uncertainty sources are correlated versus uncor- 
rected and checking that the assumption of perfect cor- 
relation is reasonable. 

Most of the equations utilized to calculate flow, 
mass, volume, and other necessary intermediate quanti- 
ties for the PVTt standard have been discussed in prior 
sections. In Fig. Al the information is summarized in a 
diagram that shows the measurement chain used to cal- 
culate flow. At the top of the diagram is the output, 
mass flow. At the second level are the inputs to Eq. (2), 
the quantities needed to calculate mass flow: density, 
volumes, collection time, etc. To calculate density, the 
inputs to Eq. (3) are necessary: pressure, temperature, 
compressibility, the universal gas constant, and molec- 
ular weight. The other necessary quantities and their 
basis equations are shown as well. Figure Al will serve 
as a guide for the PVTt uncertainty analysis. 

The discharge coefficient resulting from a flowmeter 
calibration will have additional uncertainties not con- 
sidered herein due to measurements associated with the 
meter under test. For instance, if the meter under test is 
a critical flow venturi, uncertainties related to the tem- 
perature and pressure measurements at the meter must 
be included in the uncertainty of the discharge coeffi- 
cient. 

The uncertainties tabulated herein are k=\, stan- 
dard, or 68 % confidence level uncertainties. At the 
conclusion of the uncertainty analysis, a coverage fac- 




Fig. Al. The chain of measurements and equations used for 
the Pra flow standard. 
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tor of 2 will be applied to give an expanded uncertain- 
ty for mass flow measurements with an approximate 
95 % confidence level. In the remainder of Sec. 6, the 
uncertainty components that contribute to the PVTt 
mass flow measurement will be traced to their most 
fundamental sources. 



6.1 Pressure 

A Ruska Model 2465-754' gas lubricated piston 
pressure gauge is used as the primary pressure standard 
to calibrate pressure transducers within the Fluid Flow 
Group. The uncertainties in a single pressure measure- 
ment made with this device are listed in Table Al. 
Uncertainties in the pressure standard can be traced to 
the effective area of the piston, piston thermal expan- 
sion, the masses, local gravity, and the measurement of 
the pressure under the bell jar covering the piston and 
masses (necessary for absolute pressure measure- 
ments). The uncertainty shown in Table Al is for a 
pressure value of 100 kPa. No buoyancy corrections are 
made to the masses since the reference pressure, and 
hence the density under the bell, are small enough that 
the buoyancy corrections (and uncertainties) are negli- 
gible («1 X 10"^ P). The uncertainty of the piston pres- 
sure gauge is 17 X 10"^ P at 100 kPa. 

6.1.1 Collection Tank Pressure 

Measurements of the collection tank initial pressure 
are made with a pair of thermocouple vacuum gauges 



(Varian Convectorr P-type) that have been calibrated by 
comparison to a reference standard in the NIST 
Pressure and Vacuum Group. The manufacturer's 
uncertainty specification for this gauge is 10 % of read- 
ing. Based on the NIST calibration results, the consis- 
tent agreement between the redundant sensors, and the 
repeatable readings of the gauges at the vacuum pump 
ultimate pressure, a standard uncertainty of 5 % of 
reading will be used. As will be seen when the compo- 
nents are combined to give the flow measurement 
uncertainty, this large value has little impact due to the 
low initial pressure in the collection tank (20 Pa). 

Pressure measurements of the full collection tank are 
made with a Paroscientific Model 740 with a full scale 
of 200 kPa. The manufacturer's uncertainty specifica- 
tion for this transducer is 0.01 % of full scale, but under 
the conditions of the present usage, the uncertainty is 
less. The uncertainties in the collection tank pressure 
measurement are listed in Table A2. They include the 
uncertainties from the piston pressure gage, the long 
term drift of the Paroscientific transducer which has 
been quantified by periodic re-calibrations, as well 
as the residuals from the best fit calibration equation 
(including hysteresis), and thermal effects. 
Uncertainties due to spatial non-uniformity of pressure 
within the tank and time response of the sensor are neg- 
ligible since the calibration procedure is to wait as 
much as 20 min for equilibration before the measure- 
ments are made. 

Figure A2 is a control chart that shows the changes 
in pressure calibration (xlO"*) vs time at a pressure of 
100 kPa for one of the pressure transducers used to 



Table Al. Uncertaimies for a 100 kPa pressure measurement made with the piston pressure gage used as the standard for pressure calibrations 



Uncertainty category 



Piston Pressure gage 



Standard uncertainty 
(i=l) 



Relative (xlO ) 



Absolute (kPa) 



Contrib 
(%) 



Comments 



Pressure value 

Piston effective area 12 

Thermal expansion 6 

Masses 1 

Local gravity 0.2 

Air density for buoyancy 

Ref P for absolute 10 

RSS 17 



100 

0.0012 

0.0006 

0.0001 

0.00002 



0.001 

0.0017 



65 

4 







31 



From NIST Pressure and Vacuum Group cal. 

Assume T imc of 0.2 K 

Mass Group calibration 

9.801011+/-0.000002 

Neg. air density relative to mass density 

Based on calibration data of the vac gauge 



Certain commercial equipment, instruments, or materials are iden- 
tified in this paper to foster understanding. Such identification does 
not imply recommendation or endorsement by the National Institute 
of Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the pur- 
pose. 
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Table A2. UncertaiMies in the collection tank pressure measurement at 1 00 kPa 



Uncertainty category 
Pressure measurement 



Standard uncertainty Contrib 

{k=\) 
Relative (xlO*) Absolute (kPa) (%) 



Comments 



Pressure value 

Piston pressure gage 17 

Drift 60 

Residuals, hysteresis, thermal effects 14 

RSS 64 



100 






0.0017 


7 


From Table Al 


0.0060 


88 


<0.01 % in 6 months, assume r 


0.0014 


5 


From cal. records, experiments 


0.0064 
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a, 
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a- 
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Fig. A2. A calibration control chart for a 200 kPa pressure transduc- 
er used to measure the collection tank pressure. 

measure collection tank pressure. Also shown are the 
k=\ uncertainty tolerance bounds (64 x 1 0~^ P from 
Table A2) and error bars that represent the A: = 1 uncer- 
tainty of the piston pressure gauge used to calibrate the 
sensor (17 X 10"^ P from Table Al). 

Temperature effects as large as 40 x \Q~^ P have 
been observed in the tank pressure sensor and care was 
taken to minimize their influence. When the tank is 
quickly filled from a pressurized cylinder during the 
volume determination process, cold gas enters the sen- 
sor, cooling it. The pressure readings asymptotically 
approach a final value as the sensor returns to room 
temperature (with a time constant of approximately 



1 hr). The temperature dependence of the sensor was 
confirmed by testing with an environmental chamber. 
Temperature effects also result in a hysteresis loop for 
the sensor calibration data that enlarges the calibration 
fit residuals. The calibration process entails increasing 
and decreasing the pressure in steps. The pressure steps 
result in heating and cooling of the pressure sensing 
elements and a hysteresis loop. Therefore the residuals 
of the pressure calibration fit include contributions due 
to thermal effects. We noticed that the thermal effects 
due to pressure changes in the transducers can be larg- 
er than the values given in Table A2. During volume 
determinations, we allowed sufficient time for the sen- 
sor to return to room temperature so that the remaining 
temperature effects were much smaller than the 
allowance for calibration drift. 



6.2 Temperature 

The temperature sensors used in the flow standard 
are traceable to the NIST Thermometry Group through 
calibrations made with a four-wire thermister transfer 
standard (Thermometries Model TS8901) and a recir- 
culating water constant temperature bath. The uncer- 
tainty of the transfer standard thermister is 1.2 mK (see 
Table A3). The drift of the transfer standard is consid- 



Table A3. Uncertainties for the Fluid Flow Group temperature transfer standard 



Uncertainty category 




Standard uncertainty 


Contrib 








{k = 


1) 




Temperature transfer standard 


Relative (xlO") 




Absolute (mK) 


(%) 


Temperature value 








3x10^ 




Thermometry Group cal 


4 






1.2 


94 


Fit residuals 


1 






0.3 


6 


Drift 















Radiation, self-heating, etc. 
RSS 



4 







1.2 






Comments 



Unc. for 274 K to 368 K 

Some years, 1/6 this size 

Less than discemable given cal unc. 

Deeply immersed in water bath 
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ered negligible based on seven annual calibrations that 
have always differed from each other by less than the 
calibration uncertainty. 

6.2.1 Collection Tank Temperature 

The measurement of the temperature of the gas in the 
collection tank has additional uncertainties that are list- 
ed and quantified in Table A4. Temperature is measured 
with YSI Model 46000 thermisters in 3 mm diameter 
stainless steel sheaths, a Keithley model 224 current 
source, a Keithley model 7001 switch system, and a 
Keithley model 2002 multimeter Uncertainty sources 
include the temperature transfer standard covered by 
Table A3, the uniformity and stability of the water bath 
used to calibrate the thermisters, and the residuals of 
the best-fit equation to the calibration data. The largest 
uncertainty component is the calibration drift between 
periodic calibrations. Radiation and stem conduction 
are negligible since the thermisters are immersed at 
least 15 cm in room temperature water Tests were con- 
ducted to measure the significance of self-heating by 
varying the thermister current while the sensor was 
held in a stable water bath and watching the resulting 
change in sensor reading. Based on this experiment, the 
current through the 5000 Q. thermisters was set to 
10 |iA which leads to self-heating of less than 1 mK. 
The PVTt bath stability and uniformity (1 mK) were 
discussed earlier, as was the issue of thermal equilibri- 
um between the water bath and the gas in the collection 
tank. The sensors are calibrated over their entire range 
of usage, so there is no uncertainty related to extrapola- 
tion of their calibration data. Uncertainty related to the 
time response of the thermisters is negligible since the 
time constant for the sensor is on the order of 1 s and 



the wait for thermal equilibrium is 30 or more times 
longer 

Small portions of the gas collection tank are not 
immersed in the water bath. They include the tubing 
connecting the outlet of the diverter valve to the collec- 
tion tank, the tubing that connects the tank to pressure 
and vacuum transducers, and the internal volume of 
these transducers. Because we assumed that the bath 
temperature represents the gas temperature, and the 
room temperature may differ from that of the bath, the 
small portion of the tank not immersed leads to a gas 
temperature uncertainty. The room temperature is 
maintained at 23.5 °C + 1 °C. The fractional error in 
mass contained in the collection tank due to a 1 K dif- 
ference between the room and the water bath is: 



5m 

m 



V^ T ' 



(A4) 



where V^ is the volume of the portion of the tank that 
is at room temperature and 5T'\s the difference between 
the room and water bath temperatures. For the large 
tank, Fgu, is 200 cm-' and the total tank volume is 677 L. 
For 5T of 1 K, the relative mass uncertainty is 1 x 10"^. 
This uncertainty will be treated as an uncertainty in the 
average gas temperature of 1 x 10"^ T or 0.3 mK in the 
following analysis. For the small tank this temperature 
uncertainty is 1 .2 mK or 4 X 10"^ T. 

By RSS of the components listed in Table A4, the 
combined uncertainty for the average temperature of 
the gas in the collection tank is 7 mK when using the 
thermisters dedicated to the flow standard. 



Table A4. Uncertaimy of average gas temperature in the collection tank with the dedicated temperature sensors 



Uncertainty category 




Standard uncertainty 


Contrib 


Comments 






{k- 


= 1) 






Tank gas average temperature 


Relative (xlO") 


Absolute (mK) 


(%) 




Temperature value 






297000 






Temperature transfer standard 


4 




1.2 


4 


From Table A3 


Cal. bath uniformity and stability 


3 




1 


2 


Based on notes made during cal 


Fit residuals 


7 




2 


9 




Drift (I, R, DMM, thermistors) 


20 




5.8 


77 


YSI spec is <10mK/10 months, assume rect. 


Radiation, stem cond., self-heating 


3 




1 


2 


Expt. varied current 


PVTt bath uniformity and stability 


3 




1 


2 


See Fig. 4 


Diflf. between water and room T's 


4 




1.2 


3 


See Eq. (A4) 


Extrapolation 













Cal covers whole operating range 


Sensor time constant 













Wait time is >10 x sensor time constant 


RSS 


22 




7 
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6.3 Compressibility, Molecular Weight, and Gas 
Constant 

The remaining contributors to gas density uncertain- 
ty are the compressibihty, the molecular weight, and 
the gas constant. Four gas cases will be considered. 
Tank volume determinations were conducted with ultra 
high purity nitrogen (99.999 % N2) and ultra high puri- 
ty argon (99.999 % Ar). During normal calibrations, 
industrial liquid nitrogen vaporized from dewars 
(99.998 % N2) and compressed, dried air are used. 

6.3.1 Gas Constant 

The universal gas constant is known with relative 
standard uncertainty better than 2 x 10"^ [12]. 

6.3.2 Compressibility 

The compressibility factor can be calculated from the 
following expression: 



Z = \ + Bp„+Cp„\ 



(A5) 



where p„ is the molar density (mol/cm^) and B and C 
are the second and third virial coefficients respectively. 
For nitrogen, the virial coefficients can be calculated 
from the correlations: 

B^ = -131.21 +0.65125r-7.636xl0-'r', (A6) 



Q =3454.2-11.35r+0.015r'. 



(A7) 



where temperature is in K and B and C are in cmVmol 
and cm^/moP, respectively, and the units of the con- 
stants have been surpressed. For dry air the correlations 
for the virial coefficients are: 



These expressions are the result of least squares best 
fitting to outputs from a database of the property meas- 
urements by numerous experimenters [13]. Equations 
(A6) through (A9) were fitted to data over the range of 
270 K to 330 K. This range allows application of the 
correlations at the conditions found in the test section 
of the flow standard. However, it should be noted that 
the conditions in the collection tank are much narrow- 
er, with pressures ranging from kPa to 101 kPa at a 
nearly constant temperature of 296.5 K. For these nar- 
rower conditions, the third virial coefficient could be 
ignored with negligible impact on the density uncer- 
tainty. 

Uncertainty estimates for experimental studies of 
compressibility are often unavailable, especially for 
older publications. Comparison of previously compiled 
compressibility values obtained by various researchers 
[14,15,16,17] shows agreement to within 10 X 10"^ in 
the 270 K to 330 K temperature range at 100 kPa. 
Perhaps more valuable is that this level of agreement is 
achieved between compressibility measurements made 
by the traditional /iFr method and by the more recent 
speed of sound techniques [18]. Based on this informa- 
tion, a relative standard uncertainty of 10 x 10"^ will be 
used for the experimental measurements of compress- 
ibility. This uncertainty is for a pressure of 100 kPa and 
it scales with density. The residuals of the equation fit- 
ting process to the experimental data are negligibly 
small (<1 X 10~* Z). The uncertainty components of the 
compressibility factor are listed in Table A5. 

The source of the dry air used in the flow standard is 
a Joy oil-free two stage reciprocating compressor and a 
Zum refrigeration dryer. The mole fraction of water in 
this air source is 0.0010 ± 0.0005. This uncertainty in 
composition leads to uncertainty in the compressibility 
of only 1 X 10"^ Z. The uncertainty in compressibility 
due to impurities is smaller for the other gases. 



S^ =-137.06 +0.66785 r-7.833xlO"^r', (A8) 
CAi^= 3528.9 -12.40r+0.016r', (A9) 



6.3.3 Molecular Weight 



The departure of the molecular weight of ultra high 
purity nitrogen, industrial liquid nitrogen, and ultra 

Table A5. Uncertainty of the compressibility for nitrogen, argon, and dry air 



Uncertainty category 
Compressibility (Z) 



Standard uncertainty 
(i=l) 



Relative (xlO ) 



Absolute 



Contrib 
(%) 



Comments 



Compressibility value 
Experimental data 
Impurity effects on Z 
Residuals of best fit 


10 

1 




RSS 


10 



1 

1 X 10" 
Ix 10^ 

1 X 10" 



99 

1 





Based on analysis of Uterature 



<1 X 10 
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high purity argon from the molecular weight of the pure 
substance was examined using the impurity specifica- 
tions of the gas manufacturer This analysis results in 
molecular weight relative standard uncertainties less 
than 1 X 10"*, but 1 x 10"* will be assumed for the 
molecular weight of nitrogen (28.01348 g/mol) and 
argon (39.94800 g/mol). For the dry air described 
above, the molecular weight is 28.9537 g/mol and its 
relative standard uncertainty is 190 x 10"* due to the 
variability of water content. 



6.4 Density 

Now that the sub-components have been quantified, 
the uncertainty of density measurements made in the 
collection tank with nitrogen and argon (for volume 
determinations and for flow measurements) and with 
dry air (for flow measurements) can be calculated and 
they are presented in Tables A6 and A7. For the pure 
gases, the relative standard uncertainty is 68 x 10"* and 
the primary contributor is the pressure measurement. 
For air, the relative standard uncertainty is 208 X 10"* 
and the largest contribution is due to water content vari- 
ations. The uncertainties of the density of ambient air 
(needed for buoyancy corrections) and of distilled 
water are considered in the section pertaining to tank 
volume determinations. 



6.5 Collection Time 

As explained in prior sections, the collection time is 
an approximate time measured by timers that is then 
corrected via analysis of records of the inventory pres- 
sure data and trigger voltages to minimize inventory 
mass and improve uncertainty cancellation. The base 
time is measured redundantly with two Hewlett- 
Packard Model 53131A counters. The counter calibra- 
tion and usage leads to less than 0.01 ms uncertainty. 
Due to the 3000 Hz recording frequency of the inven- 
tory pressure and trigger data, the actual rise of the trig- 
ger voltage can be any time within a 0.33 ms window. 
Assuming a rectangular distribution, the post-process- 
ing corrections will have a standard uncertainty of 0.19 
ms. This uncertainty applies to both the start and stop 
times. The time uncertainties of the two pressure meas- 
urements used in the mass cancellation procedure are 
negligible since the times are found by interpolation of 
the data records and are much smaller than 0.33 ms. 
The combined collection time uncertainty is 0.3 ms 
(see Table A8). 



6.6 Volume of the 677 L Collection Tank 
(Gravimetric Method) 

The uncertainty of the determination of the collec- 
tion tank volume by the previously described gravimet- 



Table A6. The uncertainty of collection tank gas density for nitrogen and argon 



Uncertainty category 

Collection tank density (Nj & Ar) 



Standard uncertainty 
(^=1) 



Relative (xlO°) 



Absolute 



Contrib 
(%) 



Comments 



Pressure 


64 


Temperature 


22 


Compressibility 


10 


Molecular weight (purity) 


1 


Gas constant 


2 


RSS 


68 



6.45xlO"-'kPa 




88 


From Table A2 


6.47 mK 




10 


From Table A4 


1 X 10"^ 




2 


From Table A5 


2.80 X 10"^ g/mol 









1.41xlO"^(cmV(s 


^K) 







7.82xl0"*g/cm^ 









Table A7. The uncertainty of collection tank gas density for dry air from the NIST Fluid Flow Group small compressor and dryer system 



Uncertainty category 




Standard uncertainty 


Contrib 


Comments 








{k = 


= 1) 








Collection tank density (air) 


Relative (xlO") 






Absolute 


(%) 




Pressure 


64 








6.45 X 10"^kPa 


9 


From Table A2 


Temperature 


22 








6.47 mK 


1 


From Table A4 


Compressibility 


10 








1 X 10"^ 





From Table A5 


Molecular weight (purity) 


190 








5.5xlO"Vmol 


89 




Gas constant 


2 








1.41xlO"^(cmV(s^K) 







RSS 


208 








2.38xl0"^g/cm^ 
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Table A8. Collection time uncertainties 



Uncertainty category 

Collection time 

Time value 

Timer cal and usage 

Inventory correction (start) 

Inventory correction (stop) 

RSS 





Standard uncertainty 


Contrib 




(k=l) 




Relat 


ive (xlO^) Absolute (ms) 


(%) 




100 000 







0.01 





2 


0.19 


50 


2 


0.19 


50 


3 


0.3 





Comments 



Base time uncertainty 

3000 Hz, rectangular distribution 

3000 Hz, rectangular distribution 



ric method is traceable to the uncertainty of the mass 
and density measurements made during the process, 
which are in turn dependent on the quantities shown in 
Fig. Al. The uncertainty of the density measurements 
for the pure nitrogen and argon gases used for the vol- 
ume measurements was given in Table A6. The scale 
used was a Mettler-Toledo model PR 10003, which has 
a 10 kg capacity. The uncertainty of the mass measure- 
ments of the high pressure cylinder before and after dis- 
charge is dependent on the buoyancy corrections 
(ambient air density and cylinder external volume), the 
reference masses used with the mass comparator, and 
the performance of the mass comparator The measure- 
ment of the external volume of the high pressure cylin- 
der via the Archimedes principle was described in an 
earlier section. The uncertainty of this measurement is 
traceable to the density of distilled water, ambient air, 
reference masses, and the performance of the mass 
measuring systems. 

The uncertainty of the density of ambient air during 
the course of the various weighing procedures is given 
in Table A9. The pressure, temperature, and relative 
humidity uncertainties account for instrument calibra- 
tion uncertainties as well as variations in the room con- 
ditions during the time needed to make a mass meas- 
urement. The relative combined uncertainty of the 
ambient air density is about 500 x 10"*, with the largest 
contribution being from temperature. 



The uncertainties related to the determination of the 
external volume of the high pressure cylinder are listed 
in Table AlO. The 1x10"* relative standard uncertain- 
ty for the apparent mass in air is based on uncertainty 
calculations for the true mass of the weighed cylinder 
(discussed later). The uncertainty of the apparent mass 
in water is based on the recorded temperatures of water 
which made the cylinder barely sink or barely float. 
Uncertainty of the water temperature for neutral buoy- 
ancy of 0.2 K leads to a value for the uncertainty of the 
apparent mass in water. This uncertainty is not reported 
in a relative form since the apparent mass in water is 
zero and the relative uncertainty is therefore undefined. 
The RSS is calculated using the results in the third col- 
umn and their sensitivity coefficient. The water density 
uncertainty includes thermometer uncertainties, esti- 
mates of non-uniformity of the water temperature, and 
uncertainty of the water density correlation obtained 
from the literature. As previously stated, different val- 
ues of external volume were used for the frill and empty 
cylinder due to pressure dilation, but an uncertainty 
related to pressure and temperature effects is included 
here. The collection tank volume is not very sensitive to 
this external tank volume, hence larger uncertainty val- 
ues than 51 X 10"* V^^ could certainly be tolerated. 

The relative standard uncertainties of the cylinder 
mass measurement (Table All) are quite small 
(1 X 10"*) with the major components being the refer- 



Table A9. Uncertainty of the density of ambient air 



Uncertainty category 
Ambient air density 




Standard uncertainty 
(i=l) 
Relative (xlO*) Absolute 




Contrib 
(%) 


Comments 


Pressure 
Temperature 
Relative humidity 
Equation of state (Z, M, 
RSS 


R) 


129 

336 

189394 

200 

513 








0.013 kPa 
100 mK 
2.5 % 
0.0002 
6.02 X 10"^ 


g/cm^ 


6 

45 
15 
33 


Cals and change during meas 
Cals and change during meas 
Cals and change during meas 
Based on literature references 
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Table AlO. Uncertainty of the external volume of the high pressure cylinder via Archimedes principle 



Uncertainty category 
External tank volume 


Standard uncertainty 
(i=l) 
Relative (xlO*) Absolute 


Contrib 
(%) 


Comments 


Apparent mass in air 

Apparent mass in water 

Density of air 

Density of water 

Expansion 

Std deviation of repeated meas. 

RSS 


1 

513 

20 

20 

30 

51 




4.61 X 10"^ g 
0.14 g 

6.01 X 10"^ gW 
2x10"^ g/cm^ 
9.38 X 10"^ cm^ 
0.14 cm^ 
0.24 cm^ 




34 



15 

15 

35 


Based on water r uncertainties 
Correlation and meas unc. 
Correlation and meas unc. 
T and P effects 


Table All. Uncertainty of the high- 


■pressure 


cyUnder mass measurement 






Uncertainty category 

Final mass of high P cylinder 


Standard uncertainty 
(^=1) 
Relative (xlO^) Absolute 


Contrib 
(%) 


Comments 


Reference masses 

Room air density 

Reference mass density 

External cylinder volume 

Std deviation of repeated meas. 

RSS 


0.5 

513 



51 

1 

1 




1.91 xlO"^g 

6.01xl0"''g/cm^ 
g/cm^ 
0.24 cm^ 
0.002 g 
0.00376 g 


26 

45 



1 

28 


From NIST Mass Group cal report 
From Table 10 



ence masses, the room air density (for buoyancy cor- 
rections) and the performance of the comparator 
(repeatability of the five measurements made with the 
comparator, a type A uncertainty). The complete mass 
measurement process was repeated several times for 
the full and empty cylinder conditions to assess the 
repeatability of the process, and these repeated mass 



values never differed by more than 1 x 10"' 



This 



repetition was undertaken since a previously used 
cylinder showed changes over time, probably due to 
absorption of water from room humidity variations. 
The uncertainty of the density of the reference masses 
(7.8 g/cm^) is negligible since the sensitivity of mass 
measurements to this component is extremely small. 
While some of the mass measurement uncertainties for 

Table A12. Uncertainty of the 677 L collection tank volume 



the full and empty cylinder are correlated and this could 
be used to reduce the uncertainty of the mass delivered 
to the collection tank, this benefit was not utilized since 
the improvement was not significant. A table of uncer- 
tainty components of the empty cylinder mass measure- 
ment is not presented because it is so similar to Table 
All. 

Finally, the uncertainty of the collection tank volume 
is presented in Table A 12. The relative value of the ini- 
tial collection tank density uncertainty is quite large, 
but the sensitivity of the result to this quantity is very 
small. Recall that uncertainty due to the effects of room 
temperature variations on the portion of the collection 
tank not submerged has already been incorporated as a 
temperature uncertainty. The effects of pressure 



Uncertainty category 




Standard uncertainty 




Contrib 


Comments 






(k = 


^1) 










Collection tank volume 


Relative (xlO°) 




Absolute 




(%) 




Initial mass of high P cylinder 


1 






0.00385 g 




2 


See Table All 


Final mass of high P cylinder 


1 






0.00376 g 




2 


See Table All 


Initial collection tank gas density 


lOOOOC 


) 




1.14x10"' 


g/cm^ 







Final collection tank gas density 


68 






4.25 X 10"* 


g/cm^ 


92 


See Table A6 


Expansion due to P and T 









Ocm^ 







<1 X 10"* 


Extra volume uncertainty 


1 






0.50 cm^ 







Related to Uq transfer meas. 


Std deviation of repeated meas. 


16 






10.85 cm^ 




5 


6 measurements, 2 gases 


RSS 


71 






48.44 cm^ 
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changes from vacuum to 100 kPa on the tank volume 
have been considered analytically and found to be neg- 
ligible. A small volume (1 cm') of connecting tubing 
and valve was necessary to introduce gas from the pres- 
surized cylinder to the collection tank. This volume was 
measured with alcohol and a graduated syringe and the 
uncertainty of this volume correction was estimated to 
be 0.5 cm''. By far the largest contribution is from the 
final gas density, and this in turn is nearly completely 
traceable to the pressure measurement. 

Each gravimetric volume determination required one 
day to complete due to the time required to achieve ulti- 
mate vacuum in the tank (1 Pa), the time for the pres- 
sure transducer to reach thermal equilibrium after fill- 
ing, and the time to take multiple cylinder mass meas- 
urements separated by an hour each time. The volume 
measurements were done six times, four times with 
nitrogen and twice with argon. The standard deviation 
of these six determinations was 16 x 10"* Vj, quite good 
when one considers that two different gases were used 
and the uncertainty of their equations of state. The com- 
bined uncertainty of the 677 L collection tank volume 
is 71 X 10"* Vj. The 677 L volume was used as a refer- 
ence to determine the remaining three unknown vol- 
umes via the volume expansion method. 



6.7 Volume of the 34 L Collection Tank (Volume 
Expansion Method) 



Significant uncertainties of the volume expansion 
method are related to density and the measurement of 
pressure. Many components of the previously given 
pressure uncertainty are correlated for a pressure 
change measurement made with the same pressure 
transducer over a short time period (about 1 h). 
Uncertainties of the piston pressure gauge and the sen- 
sor drift can be considered correlated for the short time 
period and small pressure changes involved (5 kPa to 8 
kPa out of 200 kPa frill scale). 

The uncorrected uncertainties in the measurement of 
pressure change are due to sensor non-linearity, resolu- 
tion, hysteresis, and thermal effects. In order to quanti- 
fy these uncertainties for the volume expansion 
method, the piston pressure gauge was used to repeat- 
edly provide a reference step change of 6 kPa to the 
pressure transducer The pressure change measured by 
the pressure transducer was compared to the pressure 
change calculated from the pressure standard. Based on 
these experiments, the uncertainty of the pressure dif- 
ference measurement during the 34 L tank volume 
determination is 0.3 Pa 

The 34 L tank volume uncertainty is summarized in 
Table A13. The largest uncertainty contributions are 
from the 677 L volume uncertainty, the density change 
in the 34 L volume, and the standard deviation of the 
repeated volume measurements. The uncertainty of the 
density change in the 677 L volume is dominated by 
uncertainty in the measurement of the pressure change. 



As explained in a prior section, the volume expan- 
sion method is performed by pressurizing a known vol- 
ume with pure gas and opening a valve to expand the 
gas into a previously evacuated unknown volume. The 
density changes and known volume are used to calcu- 
late the unknown volume via Eq. (7). In this case the 
677 L volume was used to determine the large invento- 
ry volume and the 34 L collection tank volume. 
Subsequently, the 34 L collection tank volume was 
used to determine the small inventory volume. 



6.8 Inventory Volume 

The mass change in the inventory volume is negligi- 
ble since we used the mass cancellation procedure. 
However, since there are imperfections in the proce- 
dure, uncertainty components related to the inventory 
mass change cannot be neglected. Fortunately, the most 
significant of these uncertainty components (due to 
sensor time constants) are correlated between the start 



Table A13. Uncertainty of the 34 L tank volume determined by the volume expansion method 



Uncertainty category 




Standard uncertainty 




Contrib 


Comments 








(k- 


= 1) 








Collection tank volume 


Relative (x 


10") 


Absolute 




(%) 




677 L volume 


71 






48.4 cm' 




38 


From gravimetric determination 


Large inventory volume 


1788 






1.298 cm' 







From volume expansion method 


677 L density change 


40 






3.53 X 10"' 


g/cm' 


12 


Dominated by pressure 


34 L density change 


65 






-1.12x10" 


■'g/cm' 


32 


Dominated by pressure 


Extra volume uncertainty 


15 






0.5 cm' 




2 


Measured by liq transfer 


Std deviation of repeated meas. 


46 






1.58 cm' 




16 


1 3 volume determinations 


RSS 


116 






3.95 cm' 
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and stop diversions for the methods of operation used 
in this flow standard. Other correlated inventory vol- 
ume uncertainties include the pressure and temperature 
sensor calibrations and the differences between sensed 
and stagnation values of pressure and temperature. For 
instance, the inventory temperature is measured incor- 
rectly low by the same amount at both the start and stop 
conditions due to slow sensor response time, so cancel- 
lation of temperature uncertainty occurs. More precise- 
ly, the uncertainty of the mass change within the inven- 
tory volume caused by the correlated pressure and tem- 
perature uncertainties can be expressed as: 



u (Aotj ) = 



PI 



V,M 



(T,y 



ZR 



uiT^)- 



±ru(lf)-^u(I>) 






■u(Tl) 



(AlO) 



wihjL^. ill litis uv(iiiiliuii, i/{t ^ i^^t), iJiiHvl «( irn^) ll15,^ tU^ 
uncertainties of the inventory pressure, the inventory 
temperature, and the inventory mass change during the 
collection, respectively. Note that if the uncertainties 
and the initial and final conditions are equal (i.e., 

M(r/) = M(r/), m(^) = «(^), t:=t^, and i?=^), 

then the terms within parentheses cancel, and the flow 
uncertainty related to the inventory volume is zero. 
Equation (AlO) demonstrates the benefit of matching 
the initial and final inventory conditions to optimize the 
cancellation of correlated uncertainties. 

Not all of the measurement uncertainties of the 
inventory volume are correlated. For perfect inventory 
mass cancellation, the pressure and temperature meas- 
ured at specific locations in the inventory volume must 
exhibit perfect correlation with the pressures and tem- 
peratures throughout the inventory volume. In this way, 
the sensor readings "represenf ' the conditions through- 
out the volume and when the readings at the specific 
locations match, the conditions throughout the invento- 
ry volume are matched. Unfortunately, this representa- 
tive relationship may not exist and there may be incon- 
sistencies between the pressure and temperature fields 
between the start and stop diversions. These inconsis- 
tencies may originate Irom a change in the inventory 
wall temperature or from differences in the flow paths 
between the start and stop diversions. The spatial 
inconsistencies are uncorrected and their magnitude is 
likely a function of the mass flow. 

Another source of inventory uncertainty, alluded to 
previously, is due to imperfection in matching the stop 
diversion pressure to the start diversion pressure. 



Recall that the inventory pressure is recorded while the 
bypass valve is open (nominally the barometric pres- 
sure), and that the tank filling is stopped when the pres- 
sure reaches this same pressure. In this way, the pres- 
sure at the beginning of the dead-end time transients is 
nearly equal for both diversions and the symmetry of 
the transients is improved. At high flows, it becomes 
more difficult to match these initial pressures in the 
high speed data records and a "trigger pressure differ- 
ence" that increases with increasing flow occurs. The 
size of the trigger pressure difference can be reduced by 
using faster sensors and a fast data acquisition and 
diverter valvc vontrul syDtcm. An example of the trig- 
ger pressure difference can be seen in Fig. 7, where at 
times less than zero, the pressure traces differ by about 
1 .6 kPa. The trigger pressure difference, coupled with 
the "historical" nature of the inventory pressure meas- 
urements due to the sensor time constant, is another 
reason that matching of the pressure sensor readings 
l.-^ rtul HxSv-vuri \y (jju J I . imia-ln'n^ 1 1 llu iitliiw I cit- 
ditions in the inventory volume. Hence the trigger pres- 
sure difference is another source of inventory uncer- 
tainty that scales with the flow. 

It is difficult to assess the magnitude of the uncorre- 
lated uncertainties of pressure and temperature between 
the start and stop diversions in the inventory volume. 
Our strategy is to estimate the uncorrelated inventory 
uncertainties and then perform experiments (see Sec. 4) 
to confirm that the estimates are reasonable. We 
assumed that the uncorrelated inventory uncertainties 
scale with the flow and are essentially zero at the min- 
imum flow of each tank. We will assume values of 3 
kPa and 9 K (about 3 % of the nominal values) for the 
maximum flow of each tank. 

Uncertainties related to the fast measurement of 
pressure with a 700 kPa full scale Heise Model HPO 
sensor in the inventory volume are listed in Table A 14, 
separated into correlated and uncorrelated components. 

The uncertainties of the measurement of temperature 
in the inventory volume are listed in Table A15, again 
divided into correlated and uncorrelated components. 
The uncertainties include calibration, time response, 
sensed versus stagnation issues, sensor repeatability, 
and spatial non-uniformity or inconsistency between 
start and stop diversions. 

The uncertainties for the large and small inventory 
volumes determined by the volume expansion method 
are presented in Tables A16 and A17. Although the rel- 
ative uncertainties of these volumes are large, the sen- 
sitivity of the flow measurement to these volumes is 
small due to the inventory mass cancellation scheme 
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Table A14. Uncertainties in the inventory pressure measurement, correlated and uncorrected, for maximum flow conditions 



Uncertainty category 

Inventory pressure measurement 


Relative (> 


Standard uncertainty 

(^=1) 
clO*) Absolute (kPa) 


Contrib 
(%) 


Comments 


Pressure value 
Sensor calibration 
Sensor time response 
Sensed vs stagnation 
RSS (correlated) 


3000 
1350000 
710 
1350004 








100 

0.30 

135 

0.071 

135 




100 




From inv. model at max flow 
For max. flow 


Sensor repeatability 
Spatial inconsistency 
RSS (uncorrected) 


3.00E+02 
3.00E+04 
30001 








0.03 

3 

3.0 



100 


From calibration data 
Estimated 



Table A15. Inventory temperature uncertainties, correlated and uncorrelated, for the maximum flow for each tank 



Uncertainty category 
Inventory thermocouple 


Standard uncertainty 
{k=\) 
Relative (xlO^) Absolute (mK) 


Contrib 
(%) 


Comments 


Temperature value 
Sensor calibration 


337 








297000 
100 







Sensor time response 
Sensed vs stagnation 
RSS (correlated) 


168350 

50 

168351 








50000 

15 

50000 


97 



From inventory model 
For max. flow 


Sensor repeatability 
Spatial inconsistency 
RSS (uncorrelated) 


286 

30303 

30304 








85 

9000 

9000 



100 


From calibration data 
Estimated 



Table A16. Uncertainty of the 1000 cm (large) inventory volume 



Uncertainty category 




Standard uncertainty 




Contrib 


Comments 








(k- 


= 1) 










Collection tank volume 


Relative (x 


10") 






Absolute 




(%) 




677 L volume 


71 








48.4 cm^ 









677 L density change 


1039 








3.41 X 10" 


'gW 


54 




Inv. vol. density change 


120 








2.81 X 10" 


' gW 


1 




Extra volume uncertainty 


526 








0.5 cm^ 




14 




Std. deviation of repeated meas. 


799 








0.76 cm^ 




32 




RSS 


1419 








1.35 cm^ 









Table A17. Uncertainty of the 75 cm (small) inventory volume 



Uncertainty category 

Collection tank volume 

34 L volume 

34 L density change 

Inv. vol. density change 

Extra volume uncertainty 

Std. deviation of repeated meas. 

RSS 



Relativ 


Standard uncertainty 

{k=\) 

e (xlO*) Absolute 




116 
811 
120 








4.0 cm^ 
3.41 X 10"' 
2.31 X 10"^ 


g/cm^ 
g/cm^ 


1346 
511 








0.1 cm^ 
0.038 cm^ 




1661 








0.12 cm^ 





Contrib 

(%) 



24 

1 

66 

9 



Comments 
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and the relatively small size of Fj. For the very small 
pressure changes used for the inventory volume deter- 
minations, the ability of the pressure transducers to 
measure the pressure change becomes a large uncer- 
tainty issue. The standard deviation of repeated meas- 
urements is large for these volume determinations due 
to the uncertainty of the small pressure change meas- 
urements. The ratio of the known tank volume to the 
unknown inventory volume is about 500 to 1, hence the 
pressure change is only about 0.3 kPa for an initial 
pressure of 150 kPa. The uncertainty of the pressure 
change measurement (0.3 Pa) is based on the previous- 
ly described experiments using the pressure standard. 
The instrument resolution (1.3 Pa) is also a concern for 
such small pressure changes measured with a 200 kPa 
full scale transducer Although the noise in the pressure 
data is greater than 1.3 Pa, averaging the pressure 
measurements over 30 s reduces the resolution to an 
acceptable level. For the inventory volumes, the "extra 
volume" includes tubing and valves used to introduce 
gas to the system and volume changes caused by the 
actuation of the diverter valves. The extra volume cor- 
rections are based on liquid transfer into or out of the 
volume by pipette, and /or by dimensional calculations. 
Tables A18 and A19 show the uncertainty in the 
mass change in the inventory volume for the largest 
flows in the 34 L and 677 L systems. The relative 
uncertainty of Awj is undefined since the mass change 
is zero. 



The uncertainties of gas density, volume, collection 
time, and inventory mass change provided in this 
Appendix have been combined in Sec. 3 to give the 
uncertainty in mass flow. 



Table A18. Uncertainty of the inventory mass change for the 677 L system 



Uncertainty category 




Standard uncertainty 


Contrib 


Comments 








{k = 


= 1) 








Inventory mass change 


Relative (xlO") 






Absolute 


(%) 




Initial density 


42643 








6.29 X 10"^ gW 


50 


3 kPa and 9 K uncorr. inv. unc. 


Final density 


42643 








6.24 X 10"^ g/cm^ 


50 


3 kPa and 9 K uncorr. inv. unc. 


Volume 


1419 








1.35 cm^ 





See Table Al 6 


RSS 










8.42 X 10"^ g 







Table A19. Uncertainty of the inventory mass change for the 34 L system 



Uncertainty category 




Standard uncertainty 


Contrib 


Comments 








(k = 


= 1) 








Inventory mass change 


Relative (xlO") 






Absolute 


(%) 




Initial density 


42643 








6.29 X 10"^ g/cm^ 


50 


3 kPa and 9 K uncorr. inv. unc. 


Final density 


42643 








6.24 X 10"^ gW 


50 


3 kPa and 9 K uncorr. inv. unc. 


Volume 


1661 








0.123 cm^ 





See Table A17 


RSS 










6.58 X 10"^ g 
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